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Transport properties of Manganite by double exchange mechanism are considered. The system 
is modeled by a simplified double exchange model, i.e. ,the Hund interaction between the spins of 
itinerant electrons and local spins is simplified to the Ising type. The transport quantities as the 
electronic conductivity, thermal conductivity, and the thermal power are calculated by the dynamical 
mean-field theory. It is found that the transport quantities exhibits clearly the ferromagnatic phase 
transition. A comparison with experiments is also presented. 
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I. INTRODUCTION 



The interest phenomena in the family of doped 
manganese oxides Ti-xDxMnO^ has been recently 
renewed^"**. As doping x and temperature T are varied, 
these manganite show a rich variety of phases'^. Par- 
ticulary interesting is the doping region 0.1 < x < 0.3, 
where the compounds undergo a transition from either 
insulating or very high resistance metalic, paramagnetic 
(PM) phase at high temperature to a ferromagnetic 
(FM) phase at low temperature. Near the transition, 
the resistivity of the compounds changes by a orders 
of magnitude. The application of a strong magnetic 
field substantally reduces this effect, thus giving rise 
to a very large negative magnetoresistance^. Although 
the physical mechanism, responsible for this behaviour 
has been recently the subject of much discussion and 
controversy. The double-exchange (DE) mechanism^ 
still provides a well established sterting point. The DE 
model was proposed by Zenner— who considered the 
explicit movement of electrons schematically written 
as Mnl+02^.,:iiMn^+ Mn'^+0i|,3|Mn^+ where 1,2 
and 3 label electrons that belong either to the oxygen 
between manganese or to the eg level of the Mn ions. 
In this process, there are two simultaneous motions 
involving electron moving from the oxygen to the Mn^^ 
ions and other electron from the Mii?^ to the oxygen^. 
In the DE process the motion of the itinerant electron 
favors the ferromagnetic ordering of the local spins and, 
vive verse, the presence of ferromagnetic order facilitates 
the motion of the itinerant electron. Hence, only the 
z-component part of the Hund interaction between the 
local spins and spins of the itinerant electron plays 
an essential role in the DE. In this paper we study 
the transport properties by the simplified DE (SDE) 
where only z-component part of Hund interaction is 
incorporated. The transport quanities such as the dc- 
conductivity, thermal conductivity and thermal power 
are calculated by Dynamical mean- field theory (DMFT). 
The DMFT has been extensively used for inverstigating 
strongly correlated electron systems. It is based on the 
fact that the self-energy depends only on frequency in 
the infinite dimension limit. Using the DMFT these 



transport quanities can be expressed via the spectral 
function. 

We find that the SDE captures main features of the 
transport properties of manganites. These results indi- 
cate that the DE process in manganites can be studied 
by the SDE model which is much simpler than the 
full version of the DE model. This provides a starting 
point tomard to complex including various variations to 
the DE mechanism such as the randomness, charge or 
orbital ordering. 

The poster is organized as follows. In section II we 
present the SDE model and formulas for dc- conductivity, 
thermal conductivity and thermal power. In section III 
we provide the application DMFT in SDE model. Next 
section we present numerical results for the thermal 
transport illustrating the different contributions of 
Hund coupling and concentration of itinerant electrons. 
Conclusions are presented in section IV. 



II. TRANSPORT COEFFICENTS IN SDE 
MODEL 

The Hamiltonian of the SDE model is described as 
follows 



H = - 



(ij) 



(1) 



where c]^ {c„) is the creation (annihilation) operator for 
an electron itinerant at site i with spin a. The first term 
includes the hopping only between the nearest neighbour 
sites ion, tij is the hopping integral and is scaled with 
the spatial dimension d and to have a finite result in the 
limit d — > oo.^ 
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and we take = 1 as the unit of energy. In the limit 
d -~^ oo, the bare density of states of the itinerant 
electron becomes p{e) = -^e"*^ for a hypercubic lattice. 
fi is chemical potential. The last term is the Hund 
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coupling between the spins of itinerant electron and 
local spin S^. In this model, only z component (Ising 
type) is concerned. Transport coefficient are calculated 
with in a Kubo-Greenword formalism, in which the 
dc-conductivity cr, thermal power S and the thermal 
conductivity k can be determined from relavant correla- 
tion function of the current operatosiS.. We define three 
transport coefficients as L^^, L^^ = L^^ and L^^, we have 
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(2) 
(3) 
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where the transport coefficients are found from the 
analytic continuation of the relevant "polarization 
operator" at zero frequency^", that mean 
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where L*^ {v) can be calculated from the correlation 
functions 
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after replace ii'n hy + i6 {S ^ 0+), here j and jq 
are particle-current operator and heat-current operator, 
respectively. The particle-current operator is defined by 
the commutator of the Hamiltonian with the polarization 
operator R^nii-, with our model we have. 



(9) 



where the velocity operator is Wq = '\/qe{q) and cj^^. 



is Fourier transform of the cl^: c^^ — X^o 
and the energy-current is defined by the commutator of 
the Hamiltonian with the energy polarization operator 
Ri/ii (where H = J^i^i)' with this defined a heat- 
current operator is 



q.q'o- 



JHCrS{q ~ q')(Vq' + Vq)c]j^Cq^ 



(10) 



where ^(q - q') = ^ e^^^i-i')^' The important 

relations between the heat-current and the particle- 
current operator are described fully bellow. 
Substubting j and jg from © and (p^J into ©0 and 
(|H1), the transport coefficients are easy calculated in the 
infinite dimensional hypercubic. 



L^^=TY,j dep{e) J du;{^ 
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and 



L''^Tj2l dep{e) J dcj{- 



)Alie,co)u;' (13) 



Above results lfTT|l . (fT^ and (fT!^ appears Acr{e,uj) as 
spectrum function of its Green function 



Acy{e,uj) = ImGa{t,uj) 

IT 



(14) 



That ideal lead us to calculate Green function for each 
electron with spin a. In this poster. Green function were 
found for the SDE by DMF theory approximation. 



III. APPLICATION DMF THEORY IN SDE 
MODEL 

We solve the SDE model (1) by the DMFT. The 
DMFT is based on the infinite dimension limit. In the 
infinite dimension limit the self-energy is pure local and 
does not depend on momentum. The Green function of 
the itinerant electrons with spin a satisfies the Dyson 
equation 



Ga(k,ia;„) 



iuon - e(k) -t- /i - T,a{iu}n) ' 
So we have local (single-site) Green function 



GLa{i^^n) = — 2_^Go-(k, iLJ„) 



(15) 
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where uin = 7rT(2n + 1) is the Matsubara frequency, 
e(k) — —2t^^cos{ka) is the dispersion of the free 
itinerant electrons on a hypercubic lattice, T,a(iujn) is 
the self-energy. The self-energy is determined by solving 
an effective single-site problem. The effective action for 
this problem is. 

Seff = - fdr f dr'^4(r)(7-i(r-T'K(r') 



+ 2Jh / dT5"^fT4(r)c,(T) (17) 

•' a 

where Qaij ~ t ) is the Green function of the effective 
medium QaiioJn) in the time representation . 
The local Green function of the effective single-site 
problem is solely determined by the partition function. 
It can be calculted by the equation. 



IV. NUMERICAL RESULTS 

We present the DMFT results for two cases: n — 
1.0 (half filling) and n — 0.5 (quarter filling) where 
n — —■^T'Y^^ ^ImGcr{ii^n) with two values of J = 2Jh 
[J — 2 and J — A). The algorithm for determining the 
Green function is as follows, (i) Begin with the self- 
energies in each spin (S| = 0.5 and S| = 0.0). [ii) Then 
(tT?)|l is used to find the local Green function, {iii) Subsi- 
tituting Eg- in {i) and G„ were calculated in {ii) to l|18|) 
we have the effective medium Qa- (iv) Put Q„ into H19|) 
Gcr arc determined, (w) From these Ga and Qa were de- 
termined in (Hi) and together with H18|l two new S(j are 
present. Go back to step (ii) and repeat the iteration 
until convergence is reached. In all our calculations, the 
relative error for the Green function of less than 10~^ 
is used to stop iteration loop. Our results of transport 
properties in SDE is expressed follow. 

Fig.l shows the temperature dependence of the re- 
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where Zef / is the partition function 
In addition to IjKil) . the local Green function GLa{i^n) 
can be considered as the Green function of a single-site 
problem with a certain effective bare Green function 
Q„{iuj.n) and with the same self-energy I]o-(iw„), so that 
we can write 



Gl<t(«W„) 



1 



Qa^{iuin) - Y^a{i(^n) 



(19) 



From (|16|)(|18|l and H19|) we have self-consistently equa- 
tions, which the self-energy and the Green function are 
determined. Within the effective single-site problem, the 
partition function becomes 



(20) 



where the trace is taken over Sz ■ This partition function 
can be calculated exactly, this similar to DMFT solving 
the FK model^ we obtain 



21) 



where to = — 1 + 1,...| are projections of 5* on z 
axit. 

Using Ea. H18|) we obtain the local Green function 



GLaiiuJn) = r-i 

~r ya - JHcmi 



(22) 



where 



Wr, 



Zeff 



(23) 
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FIG. 1: Temperature dependence of the resistivity p and mag- 
netization m for different electron fillings: (a) (n = 0.5), (b) 
n — 1.0. The line with filled circles denotes J = 2 and the 
line with diamonds denotes J = 4 

sistivity for different electron fillings n in two different 
Hund coupHng constant J. Fig. 1(a) provides a decrease 
in the resistivity when the magnetization increases in 
the case n — 0.5. Indeed, in the limit T — > 0, or 
TO 0.25, the density of states (DOS) of spin-up electron 
yl| (e, lS) — + fj, — e) is maximum at w = 0, that mean 
conduction electron-subsystem becomes a free electron 
gas of spin-up electron. Increasing temperature. Fig. 2 
shows that, DOS of spin-up electron is still maximum at 
Lu = and system in ferromagnetic-metalic phase with 
dependence of resistivity on temperature is quaratic. On 
the other hand, in high-temperature paramagnetic phase, 
TO = as DOS of spin-up electrons and spin-down elec- 
trons is coincides with each other but their values at 
uj = are not equal zero so system in paramagnetic- 
metal phase. 

The sharp decreases of resistivity at T < Tc for all J 
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FIG. 2: Density of states (DOS) Au{uj) as a function of fre- 
quency for different temperature T in (n = 0.5) with different 
J: left for J = 2, right for J = 4. Label 1 and 2 denotes the 
DOS for spin-up and spin-down electron, respectively. 
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is cause by rapid increases of the magnetization m. This 
discontinuity in the slope dp/ dT atT = Tc is consequence 
of the dynamical mean field approach. Incorporation of 
spatial spin fluctuationsi^i will smooth out the tempera- 
ture dependence of the resistivity in the vicinity of Tc, 
but this is beyond DMFT. 

Increasing J constant, the shape of resistivity vs temper- 
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FIG. 3: Density of states (DOS) Acr{u}) as a function of fre- 
quency for different temperature T in (n = 1.0) with different 
J: left for J — 2, right for J — 4. Label 1 and 2 denotes the 
DOS for spin-up and spin-down electron, respectively. 

ature has a litter modification but the Curie-temperature 
increases. This analysis is in qualitative agreement with 
the calculations in^ with disorder strengths A = and 
the experimental data on manganites. 
As concentration of itinerant electron increases to n = 1. 



FIG. 4: Temperature dependence of the thermal conductivity 
K for different electron fillings: (a) (n = 0.5), (b) n = 1.0. 
The line with filled circles denotes J = 2 and the line with 
diamonds denotes J = 4 



Fig. 1(b) shows more puzzle dependence of function p{T) 
on Hund coupling constant than n = 0.5 case. Indeed, 
when J is small (J = 2), Fig. 3 shows that, at low tem- 
perature, the DOS of spin-up electrons and spin-down 
electrons overlap each other, that mean, there is only one 
band energy so the dependence of resistivity on temper- 
ature is in the metalic phase. But when temperature in- 
creases, in paramagnetic phase, DOS have a gap at u; = 0, 
that mean system in insulator phase with the negative 
dp/dT above Tc- When J is larger J = 4, system displays 
insulating behaviour everywhere as Fig. 3(b) shows in all 
temperature, there are a gap at = 0, except just below 
the Curie point, where the rapid increases in the magneti- 
zation can cause the resistance to drop over a small tem- 
perature range before it turn around on increases again. 
This occurs in the transition from the paramagnetic in- 
sulator to the ferromagnetic-insulator phase because the 
charge gap in the ferromagnetic-insulator is smaller than 
the charge gap in the paramagnetic-insulator. Those 
sharp in Fig. 1(b) will generically be smoothed out by 
spatial fluctuations. 

Now, we examine the thermal properties of ours system 
inclusing the thermal power S{T) and the thermal con- 
ductivity k{T) which are shown in Fig. 4 and Fig. 5. 
The thermal conductivity behaves as expected with the 
behaviors of resistivity. In n = 0.5 case, system in met- 
alic phase in all temperature so k{T) increases every- 
where in temperature range, and when n = \ the ther- 
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FIG. 5: Temperature dependence of the thermal power S for 
different electron fillings: (a) {n = 0.5), (b) n — 1.0. The line 
with filled circles denotes J = 2 and the line with diamonds 
denotes J = 4 



mal conductivity has a sharp at Tc temperature and de- 
creases in insulator phase. At low temperature, thermal 
conductivity has a peak and then reach to zero when 
T ^ 0. This calculations agree with the calculated 
electrical thermal conductivity from electrical resistiv- 



ity in a LaQ,Q7{Ca, Pb)o,33Mn03 single crystal using the 
Wiedmann-Franz lawi. 

The thermal power behaves as expected with a linear de- 
creases to zero at low temperature which agrees with the 
behaviours of S{T) vs temperature for a mixed Pb — Ca 
doped sampled. At Tc temperature, the slope of the ther- 
mal power dS/dT has a discontinuity, but S(T) does not 
change sign in our model, that because in our model the 
itinerant-electron subsystem is not a Fermi liquid and 
the derivative of the chemical potential dfi/dT does not 
change sign at T = (does not agree with prediction 



V. CONCLUSIONS 

In this poster, we have considered the transport 
properties in SDK model by employing DMF theory. 
With DMF theory we had three self-consistent equa- 
tions which is easy calculated by numerical. Although 
the manganites are too complicated a system to be 
described completely by this simple model, we still 
arrive at some useful conclusions: The first in this model 
the various phases and picture of transition phase of 
system exit. The second, qualitative behaviours of S{T) 
and k(T) agree with experiment data, special at low 
temperature. So with this simplified model qualitative 
transport properties of DE are presented. But as we 
have said, above results only are qualitative conclusion. 
Nevertheless, basic transport properties of Maganite 
are presented and DMF theory cooperating with DE 
succeed in investigating Manganite compounds. 
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